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ABSTRACT 

The pre-main sequence star V1647 Ori has recently undergone two optical/near- 
infrared (OIR) outbursts that are associated with dramatic enhancements in the stellar 
accretion rate. Our intensive X-ray monitoring of this object affords the opportunity 
to investigate whether and how the intense X-ray emission is related to pre-MS accre- 
tion activity. Our analysis of all fourteen Chandra X-ray Observatory observations of 
V1647 Ori demonstrate that variations in the X-ray luminosity of V1647 Ori are corre- 
lated with similar changes in the OIR brightness of this source during both (2003-2005 
and 2008) eruptions, strongly supporting the hypothesis that accretion is the primary 
generation mechanism for the X-ray outbursts. Furthermore, the Chandra monitoring 
demonstrates that the X-ray spectral properties of the second eruption were strikingly 
similar to those of the 2003 eruption. We find that X-ray spectra obtained immedi- 
ately following the second outburst — during which V1647 Ori exhibited high X-ray 
luminosities, high hardness ratios, and strong X-ray variability — are well modeled as 
a heavily absorbed (Nh ~ 4 xlO^'^ cm~'^), single-component plasma with characteristic 
temperatures (kTx ~ 2-6 keV) that are consistently too high to be generated via accre- 
tion shocks but are in the range expected for plasma heated by magnetic reconnection 
events. We also find that the X-ray absorbing column has not changed significantly 
throughout the observing campaign. Since the OIR and X-ray changes are correlated, 
we hypothesize that these reconnection events either occur in the accretion stream con- 
necting the circumstellar disk to the star or in accretion-enhanced protostellar coronal 
activity. 

Subject headings: stars: formation — stars: individual (V1647 Ori) — stars: pre-main- 
sequence — X-rays: stars 
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Introduction 



The young stellar object (YSO) V1647 Ori was first noted to be an erupting source when it 
brightened suddenly in November 2003, illuminating a new nebula now known as McNeil's Neb- 
ula. This deeply embedded, low-mass YSO is typically considered to be an E X Lupi- type objec t 
(an "EXor") though it shares some spectral characteristics with FU Orionis (IVacca et al.1 12004 ) . 
the prototype of a similar class of erupting pre-main sequence (PMS) stars ("FUors"). EXors are 
observed to brighten irregularly at optical wav elengths up to several times per decade; these out- 
bursts persist for weeks to months (lHerbigll200ll ). FUors are YSOs that erupt less often than EXors, 
perhaps only once per century, and fade much more slowly, i.e., on timescales of years to decades 
(|Herbigi 119771 ). The general consensus is that these eruptions are the result of massive accretion 
events that occur irregularly as young protostars grow; such accr etion episodes may be the p rimary 
process through which young stars accrete most of their mass (jHartmann &: KenvonI 119961 ). The 
mechanisms underlying EXor and FUor outbursts remain very poorly understood; however, models 
generally invoke the rapid onset of disk instabilities that lead to sudden inward migration of the 
inner disk truncation r adius and dr a matic changes in the star/disk magnetic field configuration 



(e.g., IZhu et al.l ([20091); lK6nigl et al.l (1201 ih ). 



Soft X-ray (~0.3 keV) generation via shock-heating of plasma at the bases of PMS accre- 
tion streams has be en proposed for some objects from X-ray grating observations (see review by 
Giidel &: Nazell2009l ') using plasma diagnostic s based on helium-lik e ions (see review by | Porquet et al 
20 id ) , a prime example being TW Hydrae (jKastner et al.l l2002l : iBrickhouse et al.l |2010|) • On the 
other hand, studies such as the Chandra Orion Ultradeep Project (COUP) indicate that coronal ac- 
tivity, and not accretion, is likely the primary p roduction mechanism for the harder X-ray s that are 



characteristic of most pre-main sequence stars (jPreibisch et al 



2005 



Stassun et al.l 12006 ). Debate 



still lingers as to the primary X-ray generation mechanism in YSOs during optical outbursts. 

Since 2003, V1647 Ori has been obser ved to under g o two distinct outb ursts that have been 
detected in t he optical and near - infrared (IMcNeill l2003l: lltagaki et al.l l2008l) as well as at X-ra y 



wavelengths (jKastner et al.l 12004 l2006l : ICrosso et al.l l2005l : ICrossol l2006l : iHamaguchi et al.l l2010|) . 

The magnitudes and timescales for brightening and fading displayed by V1647 Ori do not match 
those of EXors. Furthermore, the duration of the initial outburst was shorter than those of FUors 



(Kastner et al 



2006 



Aspin. Beck. &: ReipurthI 120081 ). How ever, spectroscopic data do show signa- 



tures of accretion, such as strong Ha and Br7 emission (IReipurth &: Aspin 



2004 



Brittain et al 



2010l ). confirming that the eruptions of V1647 Ori do resemble those of FUors and EXors. 



V1647 Ori has been observed in X-rays before and during both of the eruptions (as well as 
after the first of the two eruptions) detected in the optical/near- infrared s ince 2003. Althoug h 
other FU or- or EXor-l i ke Y SOs have been observed in X -rays (V11 18 Ori (lAudard et al.ll2010l ): 



Cyg (jSkinner et al 



200 



EX LuD (IGrosso et al.l 120101 '): Z CMa (jStelzer et al.ll2009l ): FU Ori (ISkinner et al 



200fi[ h V1735 



V1647 Ori is the only eruptive YSO to undergo such extensive X-ray 
monitoring. During the 2003 eruption, the sudden rise in flux and subsequent decline in both the 
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optical and near-infrared correlated str ongly with a sharp i ncrea, se and then decline in the X-ray 
flux, which suggests a common origin ([Kastner et al.ll2004l . |2006| ). In this paper, we present and 
analyze X-ray spectra obtained with the Chandra X-ray Observatory (CXO) during the 2008-2009 
outburst. In ^ we describe the observations and data reduction. In S}3l we examine the trends 
and patterns seen in the X-ray emission from V1647 Ori over the past two outbursts, present 
modeling results for each of the five recent CXO observations, and compare these results to those 
obtained from previous (2003-2005) CXO and XMM-Newton observations. Finally, in Jj4]we discuss 
the implications of the results, and we argue that the measured properties of the X-ray emitting 
plasma are best understood within the framework of the accretion-generated X-ray hypothesis. 



Observations & Data Reduction 



Observations of V1647 Ori made with the Chandra X-ray Observatory (CXO) were triggered 
in 2008 September with a 20 ks observation (CXO Cycle 10, PI: D. Weintraub, ObsID 9915), after 
V1647 Ori was report ed to have uiidergo ne a new optical outburst between early 2008 January 
and late 2008 August (lltagaki et al.l 120081 ). Subsequent 20 ks observations were initiated in 2009 
January and April (ObsIDs 9916 and 9917, PI: D. Weintraub). In addition, in 2008 November, 
two observations (ObsIDs 10763 and 8585; PI: N. Calvet) were made of the NGC 2068/2071 region 
that serendipitously included V1647 Ori. Together, these five CXO pointings yield an extended 
sequence of observations of this source over a seven-month period immediately following the onset 
of the optical outburst in 2008. 

Fortuitously, a field of view that i ncludes V1647 Ori w as observed in 2002 November, a full year 
before the start of the 2003 outburst (ISimon et al.ll2004l ): subsequent CXO observations targeting 
V1647 Ori w ere obtained in 2004 Marc h and April, 2005 August and December, and 2006 May 



and August (IKastner et al 



2004 



20061). In addition, targe ted XMM-Ne wton observations were 



obtained in 2004 April (jOrosso et al.ll2005l ) and 2005 March (|Grossdl2006l ). and a targeted Suzaku 
observation was obtained in 2008 October (|Hamaguchi et ahlboiol ). The 2004-2006 observations 
cover a two-year period during which V1647 Ori more or less faded steadily from a strong YSO 
X-ray source into a very faint one, albeit with large X-ray flux variability. For direct comparison 
with our analysis of the 2008-2009 epoch datasets, we have re-reduced and analyzed all of the 
2002-2006 CXO datasets, using the same techniques and software packages that we have used for 
the latest (2008-2009) datasets. 

For all Chandra observations, the Advanced CCD Imaging Spectrometer (ACIS) was used in 
one of two imaging configurations. ACIS detectors have a pixel size of 0.49" and each ACIS-I and 
ACIS-S CCD has a field of view of 8.3' x 8.3'. CXO/ACIS has significant sensitivity over the 
energy range 0.3-10 keV, with the soft (<1 keV) X-ray sensitivity dependent on whether a back- 
illuminated or front-illuminated CCD serves as the detector (see below). The CIAO v4.1 software 
package and CALDB v4. 1.0-4. 1.4 calibration files were used to reduce the data and extract pulse- 
invariant (PI) spectra. Observation details are given in Table [TJ To compare X-ray luminosities 
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from 2008-2009 with those obtained from previous observations, source spectra were extracted from 
2.5" radius regions (making sure the aperture was positioned so as to encompass as many source 
photons as possible, even when V1647 Ori is 4' off-axis) while background spectra were extracted 
from regions near but beyond 2.5" from the target, on the same chip, using 20" outer radius 
extraction apertures. Inspection of background light curves reveals background levels to be fairly 
constant with no evidence of large fluctuations, such as flares, occurring during the observations. 
Spectra were also re-extracted from previous Chandra observations using the same apertures to 
ensure that direct comparisons of results from different observation epochs could be performed 
reliably. For observations in which the primary target was not V1647 Ori (ObsIDs 2539, 10763, 
and 8585), exposure maps were generated to investigate off-axis position effects. In all three of these 
observations, the net count and mean count rate corrections due to off-axis source positions were 
insignificant, on the order of a few percent. The resultant spectral data points of the five recent 
observations were grouped into energy bins with a minimum of five counts per bin before spectral 
modeling was done. The count rates were high enough and durations long enough for each of the 
2008-2009 observations that this bin size yielded PI spectra with good statistics. For the 2002-2006 
observations, we employed single-count-minimum binning for very low-count observations and five- 
count-minimum binning for higher-count observations so as to yield spectra suitable for modeling. 
Because no events with energies less than 0.5 keV and few events with energies greater than 8.0 
keV were detected from V1647 Ori, we limited our modeling to the 0.5-8.0 keV energy range. 

Since V1647 Ori was imaged with both front- and back-illuminated CCDs in the exposure 
sequence under analysis, we generated synthetic spectra and convolved these spectra with the in- 
strument responses to determine whether mean, broadband (0.5-8.0 keV) count rates were directly 
comparable for all observations of V1647 Ori. The simulations (Table [2]) showed that, for the 
plasma temperature regimes considered here, the back-illuminated S3 CCD is ~10% more sensitive 
to incoming flux than the front-illuminated ACIS-I CCDs. Net counts and mean count rates in 
the energy range 0.5-8.0 keV (Tabled]) were therefore adjusted downward accordingly for the two 
observations (ObsIDs 5307 and 5308) that used the S3 CCD in order to remove this sensitivity 
bias. 



Results 



3.1. Short-Term and Long-Term Variability 



We calculated the median photon energies, mean count rates, and mean hardness ratios for 
all 14 CXO observations of V1647 Ori. These results are reported in Table [H In calculating the 
hardness ratios, the hard (H) X-rays are defined as those in t he range from 2.8 to 8.0 keV and 
the soft (S) X-rays as those with energies from 0.5 to 2.8 keV (iGrosso et al.ll2005l ). The hardness 
ratio is then HR = (H-S)/(H+S), such that negative values of HR indicate softer spectra, and 
positive values indicate harder spectra. The recalculated values and error ranges of the mean count 
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rates and m edian photon energie s for the 2004-2006 observations differ only i narginany from those 
reported bv lKastner et al.l (|2006l ). With data from only three epochs in hand. iKastner et alj (|2004l ) 
reported that the X-ray flux brightened and hardened during the outburst (2002 November 14 to 
2004 March 07) and then quickly faded and softened post- outburst (2004 March 07 to 2004 March 
22). With data from additional epochs available, however, iKastner et al.l (|2006l ) reinterpreted the 
late 2004 March data as a short-term downward fluctuation, not as a quick end to the original 
outburst. They interpret the observed pattern of changes in the median photon energy as evidence 
that the X-ray spectrum hardened during outburst, remained somewhat hard for at least one year, 
and then softened after 2005 August 27 as both the optical and X-ray flux from V1647 Ori returned 
to pre-outburst levels. The results obtained here (for hardness ratio) confirm this general trend. 
This hardening and softening (the mean hardness ratios (with la uncertainties) changed from —0.34 
lb 0.31 (soft) to -1-0.46 ± 0.15 (hard) and back to —0.45 ± 0.33 (very soft) over the course of three 
years, from 2002 November through 2005 December) took place over a timescale of months, with 
V1647 Ori remaining near peak outburst X-ray luminosity and hardness for more than one year - 
with the notable exception of the 2004 March 22 data. V1647 Ori likely remained in its soft X-ray, 
low X-ray luminosity phase for almost three years, from late 2005 until its latest optical outburst in 
2008 August. The five CXO observations covering the interval from 2008 September through 2009 
April all show V1647 Ori in a hard X-ray, high X-ray luminosity state. As was the case during the 
2002-2006 observations, the source remained near peak X-ray luminosity for at least seven months 
(Fig. [m figure to be discussed in Section H]). As of 2009 December 18, V1647 Ori appeared to still 
be in an outburst phase with J-, H-, and K-band photome tric magnitudes near values measured 
at the brightest phase of outburst during 2008 September (jVenkat &: Anandaradl2009l ). Thus, it 
appears that V1647 Ori remained near peak luminosity for at least one year.^ 

X-ray light curves and hardness ratios as a function of time were extracted for each observation 
to examine the properties of any possible short-term X-ray variability. In particular, these time 
series data permit us to investigate whether the trends in X-ray flux and hardness ratio seen in 
2002-2006 were repeated during the recent outburst that began around 2008 August. The extracted 
X-ray light curves and hardness ratio time series for six of the nine observations from 2002-2006 are 
presented in Figure [TJ X-ray light curves and hardness ratio curves for the five 2008-2009 Chandra 
observations are presented in Figures [IHSl which show (from top to bottom) the overall X-ray light 
curve, the soft (0.5-2.8 keV) light curve, the hard (2.8-8.0 keV) light curve, and the hardness ratio 
time s eries. Energy ranges for the soft and hard X-ray bands follow those defined in lOrosso et al.l 
(|2005l ). and plotted uncertainties in mean hardness ratios and count rates are la. 



During all five of the most recent (2008-2009) observations, V1647 Ori experienced very pro- 
nounced, short-duration (~2-10 ks) X-ray flux variations in the form of large rises and subsequent 
declines in count rate. These fluctuations are evident at 4 ks and 16 ks in Figure [2l at 8 ks in Fig- 
ure O at 20 ks in Figure HI at ~8 ks in Figure [5l and at the beginning of the observation reported 



^As of submission of this paper, the outburst appears to be ongoing. 
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in Figure [6j The mean count rate of V1647 Ori varied by about a factor of ~2 for the first two 
epochs, while in the last three epochs, the mean count rate varied by a factor of ~6. In all five 
epochs, the bulk of the X-ray variability is contributed by the harder portions of the spectra. 

The X-ray count rate for V1647 Ori was intrinsically lower during most of the 2002-2006 
observations (the count rate was >5 ks~^ only on 2004 March 07) than in 2008-2009 (the count 
rate was >5 ks~^ during all five observations). Also, in some cases like ObsIDs 5307 and 5308 
(both from 2004), the observations were on the order of only a few kiloseconds in duration. The 
low count rates in most of the 2002-2006 data and the short exposure times in two of these 
observations make comparisons with the more recent datasets, in terms of the presence or absence 
of 2-10 ks timescale variations, fairly difficult. However, observations 5382 and 5383, obtained in 
2005 April and August, respectively - more than a year after the initial optical outburst of V1647 
Ori in 2003 November - yielded enough counts that we can compare these two light curves with 
the more recent observations with some confidence. As seen in these two light curves from 2005 
(Fig. [TJ lower-left panel, beginning at ~10 ks, and lower- middle panel, beginning at ~5 ks), the 
X-ray count rate again varied by a factor of ~4 on a timescale of 2-10 ks, a level and duration of 
short-term variability very similar to that seen in the more recent (2008-2009) observations. We 
conclude that kilosecond-timescale variations are seen in all observations of V1647 Ori in which the 
signal is strong enough and the observing window long enough for us to search for such events, and 
that the amplitude duty cycle and duration of these few-kilosecond X-ray variations were stable 
over the five-year period from 2005 to 2009. Such variability appears to be part of the normal X-ray 
pattern for this object, at least when the count rate is elevated. One therefore needs to interpret 
the overall flux levels observed during the 2004 March 07 and 2004 March 22 observations, which, 
respectively, are only 5.5 ks and 4.9 ks in duration, with caution. 

The highest measured X-ray mean count rate following the latest (2008) optical outburst of 
V1647 Ori is higher (23.6 ks~^; one month post-outburst) than the highest level (11.7 ks~^; four 
months post-outburst) observed during its 2003-2005 outburst (Table [TJ. In addition, with the 
exception of the 2004 March 07 observation, the 2008-2009 observations (a period spanning seven 
months) show V1647 Ori to be consistently brighter than it was during the 2003-2005 outburst. 
We find that in 2008-2009 the X-ray spectrum stays relatively hard, as is shown by the median 
photon energies, which in 2008-2009 are consistently about twice the energy (3.4 to 3.9 keV) that 
was observed before the first eruption (2.0 keV in 2002) and after V1647 Ori had returned to its 
quiescent state in late 2005 (1.7 keV). Also, while the mean hardness ratios were negative both 
before (2002 November) and after (2005 December) the first eruption, the mean hardness ratios in 
2008-2009 stay consistently around ~0.3 or greater. 

These results also reveal that V1647 Ori exhibited a very interesting sequence of changes of 
state, beginning in 2004 March. V1647 Ori had switched from its elevated state on 2004 March 
07 (11.7 counts ks~^, median photon energy of 3.6 keV, mean hardness ratio of -1-0.46) to what 
appears to be a quiescent state on 2004 March 22 (2.5 cts ks~^, median photon energy of 2.0 keV, 
—0.45 mean hardness ratio) but then was again in an elevated state one year later on 2005 April 
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11 (4.8 cts ks~^, median photon energy of 3.5 keV, mean hardness ratio of +0.51). Yet by 2005 
August 27, V1647 Ori was again returning to quiescence (0.9 cts ks~^, 3.0 keV median photon 
energy, +0.01 mean hardness ratio) and by 2005 December was fuhy back to quiescence. 

At least three explanations are plausible for the observed changes we have described in the 
X-ray emission from V1647 Ori that occurred between 2004 March 07 and 2004 March 22 and 
between 2004 March 22 and 2005 April 11. One: the outbursts generate both hard and soft X-rays; 
most of the time (2004 March 07 and 2005 April 11), the hard X-ray plasma is dominant in the 
X-ray spectrum as observed by CXO, but at certain times (2004 March 22) the soft X-ray plasma 
component is seen more clearly. Two: the first major outburst (2003 November through early 
2004 March) had ended, at least in the X-ray regime, after only a few months. By 2005 April, the 
source was again in an elevated state and by 2005 August, four months later, was again returning 
to quiescence. Three: the 2004 March 22 observation caught V1647 Ori during a short duration 
(2-10 ks) quiescent moment during the extended (months long) outburst phase. However, the 
XMM-Newton observation obtained on 2004 April 4 showed a low X-ray flux level at the beginning 
of the observation that is consisten t with the quiescen t level observed by CXO that increased in 



the second part of the observation (jGrosso et al.ll2005l ). Moreover, the supporting evidence from 



observations made at optical and near-infrared wavelengths for the first outburst ending in 2004 
mid-March and for a second outburst beginning later that year (explanation "Two") is weak, at 
best. Of the remaining two possibilities, we will argue in §4 that explanation "One" is more likely. 

A comparison of the count rates (Figures [IHSl top panels) and hardness ratios (Figures [2H6l 
bottom panels) during the ~2-10 ks variability periods suggests that the X-ray spectrum of V1647 
Ori may harden as it brightens and soften as the source fades during these few-kilosecond variations. 
This appearance of a correlation between the changing count rate and the changing hardness ratio 
on a few kilosecond timescale during short-term variability, however, is not consistent. In fact, 
we find a weak anti-correlation between these properties (correlation coefficient = —0.18) on these 
short timescales; however, we do find a strong, positive correlation over much longer periods of time 
between the X-ray luminosity (calculated in our spectral modeling; see Section 13. 2p and hardness 
ratio, as shown in Figure [71 In this figure, the top panels show the seven-year X-ray luminosity 
curve of V1647 Ori while the middle panels show the corresponding hardness ratios with a single, 
average measurement plotted for each observation epoch. A more direct demonstration of this 
correlation is shown in Figure [8] where we plot the change in mean hardness ratio versus the change 
in X-ray luminosity between observations of V1647 Ori. We derive a correlation coefficient of 0.44 
for these two parameters. 

As seen in the top and middle left-hand panels of Figure [71 from the first Chandra observation 
in 2002 to the first observation of V1647 Ori in 2004, the luminosity increased by ~2 orders of 
magnitude while the mean hardness ratio increased from a relatively soft value of —0.34 to a hard 
value of 0.46. From the second to the third observation, the X-ray luminosity dropped by an 
order of magnitude while the mean hardness ratio softened to a much lower value of —0.45, much 
lower than that seen in the 2002 observation of V1647 Ori when it was seen to be in quiescence. 
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When the fourth Chandra observation occurred in 2005, the luminosity had risen back to outburst 
levels again while the mean hardness ratio increased again to +0.51, the highest observed to that 
date. In subsequent Chandra observations through early 2006, the X-ray luminosity of V1647 Ori 
diminished while the hardness ratio dropped. 

In the top and middle righthand panels of Figure [3 we present the X-ray luminosity and 
mean hardness ratio versus time for most recent Chandra observations of V1647 Ori. The same 
correlation between changes in X-ray luminosity and changes in the mean hardness ratio seen in the 
2002 through 2006 observations is evident in all observations from 2008 through 2009. From 2008 
September 18 to 2008 November 27, there appears to be an anti-correlation with the luminosity 
decreasing slightly between these observations while the hardness ratio increases. We note, though, 
that in both of these epochs, the luminosity is high and the spectrum is quite hard. Since the short- 
period light curves of each observation show that there is typically X-ray variability on timescales of 
kiloseconds, with X-ray flux sometimes increasing by several factors, this apparent anti-correlation 
might be attributed to short-duration variability. 

These results suggest that we can identify two states for V1647 Ori, as seen in time-averaged 
data on a 5-30 ks timescale: 

• In its quiescent state, V1647 Ori has a low count rate (<1 ct ks^^), low median photon energy 
(<2.5 keV), and negative (soft X-ray dominated) hardness ratio. 

• In its elevated state, the count rate is at least 25-50 times greater, the median photon energy 
doubles, and the hardness ratio becomes strongly positive (hard X-ray dominated). 

The X-ray evidence shows that when V1647 Ori enters a major optical/near-infrared outburst, the 
X-ray profile switches from quiescent to elevated, and when the outburst phase ends, the X-ray 
profile switches from elevated to quiescent. 



3.2. Spectral Modeling 

3.2.1. Procedure & Results 

Modeling of the spectra employed XSPEC vl2.4. For those spectra with single-count binning, 
XSPEC was set to use the Cash-statistic instead of the statistic to assess the goodness of fit. All 
models at first assumed a thin, single-temperature plasma (APEC model) subject to absorption 
by an intervening column of hydrogen (WABS component). From the 2004 April XMM-Newton 
observations, single-component thermal plasma models of the V1647 X-ray spectrum yield a hy- 
drogen column density, chemical abundance, and plasma temper ature of Nj^ = 4.1x10^^ cm~^, Z 
= 0.8 solar, and kTy = 4.2 keV, respectively (jCrosso et al.ll2005l ). and recent Suzaku observations 
yield similar results (iHamaguchi et al.ll20ld ). Therefore, initially, N/f, Z, and kTx were set ap- 
proximately to the aforementioned values {^h = 4.0 x 10^^ cm~^, Z = 0.8 solar, and kTx = 4.0 
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keV). As our initial test indicated that Z is poorly constrained by the Chandra data, only Nj/ and 
kTx were left free to vary during the fitting procedure. For four of the five 2008-2009 observations, 
the models converged to best-fit solutions with column density and plasma temperature within a 
90% confidence interval of their initial values; however, the 2008 November 27 spectrum model was 
unable to converge to a physically-meaningful fit, so the plasma temperature was fixed at kT^ = 
4.0 keV (Table ED. 

A similar fitting method was used for the 2002-2006 observations (Table H]). Initially, the 
first five of these spectra spectra were modeled with freely varying hydrogen column densities and 
plasma temperatures, while the chemical abundance was fixed at 0.8 solar. These spectral models, 
however, were unable to constrain plasma temperatures or X-ray fluxes and luminosities, so the 
hydrogen column density was then fixed at Njy = 4.1 x lO^'^ cm~^, and the models were refit. The 
remaining four spectra, obtained when V1647 Ori was reverting to an optical/X-ray quiescent state, 
were modeled with the hydrogen column density and plasma temperature fixed at Njj = 4.1 x 10^^ 
cm~^ and kTx = 0.86 keV, respectively. Given the small number of counts in these four spectra, 
the associated error ranges for these X-ray fluxes and luminosities were obtained by multiplying the 
error ranges for the mean count rates of these observations with an appropriate energy conversion 
factor (ECF), where the ECF for each observation epoch was obtained by dividing the derived 
X-ray flux by the mean count rate. 

For the CXO observations of V1647 Ori with sufficient total counts (ObsIDs 9915 and 9917), 
we performed fits of a two-component thermal plasma model with parameters for hydrogen column 
density set to 4.0 x 10^^ cm~^, plasma temperatures set to kTx = 0.5 keV and 2.0 keV, and 
chemical abundance fixed at 0.8 solar. Visual inspection of the spectral models shows that there 
is a negligible difference between the best-fit single- and two-component models, F-test results 
suggest that there is no statistical improvement in the model fits with the addition of a second 
plasma component; i.e., the best-fit parameters for the latter model converge on values such that 
the contribution of the lower-temperature component is negligible. We conclude, therefore, that all 
of the 2008-2009 CXO data are best fit with a single-component model. 

Best-fit models for each of the five recent Chandra observations are shown in Figure M The 
overall trend of the spectral models is to con verge to fits with pa rameters similar to those of the 



best-fit single-component model reported bv iGrosso et al.l (j2005l ). Intervening hydrogen column 



densities do vary from model to model but remain in the ~ 4-6 x 10^^ cm ^ range, and plasma 
temperat ures are kTy ^ 2-6 keV, slightly higher than but still within the range of uncertainty 



found bv lKastner et al.l (|200g). 



Determining a robust model for the last three recent epochs was more challenging than for 
the first two epochs, especially for ObsID 9916 (2009 January 23). A major concern for us was 
to assess whether the flux and hardness variability seen in the light curves for these observations 
(Figs. SHS]) contributes to the difficulty of fitting a unique model to the spectral data. To assess 
the effects of variability on the robustness of our spectral analysis, we divided each of these the 
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spectra into two phases, one in which the X-ray flux was relatively lower and one in which the flux 
was elevated. Each of these spectra was then modeled independently. Though the hardness ratios 
suggest that the spectra did harden marginally as V1647 Ori went from X-ray-dim to X-ray-bright 
during these observations, the spectral modeling yields no measurable changes in hydrogen column 
density or plasma temperature. Thus, it appears that the changes in the light curves are likely due 
only to changes in the plasma emission measure. 



3.2.2. Iron Line Emission Analysis 

Our initial modeling of the V1647 spectra from 2008 and 2009 suggests the presence of line 
emission from near-neutral iron at 6.4 keV as well as from the (unresolved) helium-like iron Ka- 
line triplet at 6.64, 6.67, and 6.70 keV. The neutral iron line emission at 6.4 keV is often seen in 
accretion-powered sources and is usu ally attributed to the fl uorescence of cold gas in the presence of 



a nearby X-ray continuum emission (jTsuiimoto et al.ll2005l ). Thus, the presence of this line in some 



of the V1647 Ori spectra could be attributed to fluorescence of (neutral) circumstellar disk material 
by accretion-generated X-rays. It would not be surprising to detect this emission in the spectra of 
V1647 Ori given that the environment of this YSO does appear to contain the necessary ingredients 
for the formation of 6.4 keV emission, namely a strong, relatively hard X-ray source illuminating 
cold circumstellar material. Observations of other YSOs, including roughly a half-dozen sources in 
the COUP survey, with these environmental components have also shown this feature. 



Following iGrosso et al.l (|2005l ). we added a Gaussian component centered at 6.4 keV to the 
spectral models in order to account for the neutral line component. Models of two observations 
(2008 September 18 and 2009 April 21) appear to be well fit with the addition of a 6.4 keV line 
with equivalent widths of ~200 and ~500 eV, respectively, while the 2008 November 27-28 and 
2009 January 23 are fit well without the additional neutral iron line emission component (Fig. [TO]) . 
On the other hand, it is possible that the 6.4 keV emission is present during all of the 2008-2009 
observations from CXO but that its spectral signature is muffled by noise in the 2008 November 
and 2009 January observations. While it appears that we have detected 6.4 keV iron emission in 
the 2008 September 18 and 2009 April 21 observations, we cannot conclude deflnitively that these 
"detections" are real. We note that the best-fit 6.4 keV equivalent widths are poorly constrained 
with error ranges extending from zero to roughly twice the equivalent width values. Given this 
range of uncertainty, we cannot exclude the possibility that there is no 6.4 keV iron line. We feel 
confident, however, that these detections are real given that v i sual a nd quantitative comparisons 



of our findings with those detections found by lTsujimoto et al.l (|2005l ) are very similar. 



We have compared the intensities/appearances of this line in the various observations of V1647 
Ori by CXO, XMM-Newton, and Suzaku. The equivalent widths of this emission feature are very 
similar to those found in the spectra of the 2004 April 4 XMM observation (109 eV) and in the 
2008 October 8 Suzaku observation ('--'600 eV). We are unable to clearly determine whether there 
is a correlation between the strength and appearance of the 6.4 keV line and any of the associated 
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plasma characteristics. 



4. Discussion 



Figure [TT] suggests that the overah X-ray flux of V1647 Ori is strongly correlated with optical/near- 
infrared flux. This correlation is revealed more clearly in Figure [T2l in which we plot the X-ray 
luminosity versus the Ic-band luminosity. We interpolated Ic-band luminosities for the 2002-2005 
X-ray observation dates. We did not extrapolate I^-band luminosities for the two 2006 X-ray ob- 
servation dates due to the highly uncertain flux behavior of V1647 Ori. Of the 2008-2009 X-ray 
observations, we could only interpolate an Ic-band luminosity for ObsID 9916. With these eight in- 
terpolated luminosities, we derived a correlation coefficient between X-ray luminosity and Ic-band 
luminosity of 0.65. 

Dramatic increases in optical/near- infrared flu x for YSOs, such as FU Or i, have long been 
thought to be associated with enhanced accretion (jHartmann &: KenyonI Il996l ). In such an en- 
vironment, material is channeled through n iagnetic funnels from the co-rotation radius of the 
circumstellar disk down to the photosphere (jShu et al.l Il994l ). Hence, when the X-ray flux from 
a pre-main-sequence star or protostar is elevated and the rapid rise in X-ray emission is directly 
correlated with large-scale optical outbursts, the correlation itself strongl y suggests that accre tion 
is the mechanism responsible for generating the increase in X-ray output (jKastner et al.ll2006l ). 



One way in which accretion-generated X-ray emission could be identifled observationally would 
be through the relatively soft X-rays emitted by the plasma when it plunges onto the stellar surface 
at fre e-fall velocities and is shock heated to ternperatures of a few million Kelvin (k Ty ~0.3 
keV) (jKastner et al.ll2002l : Istelzer fc SchinittlbooilRaassenlbood : ISrickhouse et aLlboid ). We can 



speculate that a second such signature would be muc h hotter (kTy ~few to tens of keV) and harder 



X-ray s generated in magnetic reconnection events (IShibata &: Yokovamal 120021 : iBrickhouse et al 



20 id ) within the accretion streams. Whether the soft or hard X-ray generating plasma is observable 
likely depends on the observing geometry. If we have an unobscured line of sight to the footprint 
of the accretion column, our observations should be sensitive to the cooler plasma; if the accretion 
column obscures our view of the accretion-column footprint, our observations should be sensitive 
only to the hotter plasma in the accretion stream; and if the obscuration of the accretion footprint 
is partial, we might detect both plasmas. In between periods of dramatic ally enhanced accretion, 
previous studies such as COUP (jPreibisch et al.ll2005l : IStassun et al.ll200a) and the XMM- Newton 
Extended Survey of the Taurus molecular cloud (XEST) (jAudard et al.l 120071 ) , suggest that the 
X-ray signature of young stars should be that of normal coronal emission. Such emission would 
be similar to but much fainter than the hot, hard plasma seen from reconnection events in the 
accretion funnel. 



During quiescence, the X-ray flux of V1647 Ori has a low count rate, low median photon energy, 
and negative hardness ratio; in contrast, when the X-ray flux is elevated, the count rate increases 
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by a factor of 25 or greater, the median photon energy doubles, and the hardness ratio becomes 
strongly positive. Most of the V1647 Ori light curves reveal that this YSO also experiences what 
appear to be short-term (few kilosecond) variability in its X-ray flux. 



4.1. Observations of Hotter and Cooler X-ray-Generating Plasmas 

Since most of the CXO observations show this short-duration variability when the X-ray count 
rate is high, the short-term variations are likely part of the normal behavior for V1647 Ori. For 
all 2008-2009 CXO observations, the X-ray variability is almost entirely seen in hard (>2.8 keV) 
X-rays. Our model fits of the V1647 Ori spectra indicate that the X-ray signature of the plasma 
we observe from this source during optical eruptions is predominantly a bright, hot, hard, single- 
temperature plasma heated to temperatures of 2-6 keV , which is also cq r isistent with X-rays ge n- 



erated from reconnection events in an accretion funnel (|lsobe et al.ll2003l : iBrickhouse et al.ll2010l ). 



We must keep in mind, however, that the intervening hydrogen column density will play a 
significant role in our ability to detect X-ray fiux from accretion hotspots at temperatures of a few 
million Kelvin, especially when dealing with column densities as high as those modeled in V1647 
Ori observations. An accretion-footprint (shocked) plasma at a temperature of kTx ~ 0.3 keV 
subject to an intervening absorbing column density of a few times 10^^ cm~^ would have ~98% of 
its X-ray fiux extinguished. Even given the smallest hydrogen column density found via spectral 
modeling of the 2004-2005 VI 647 Ori observations when the hydrogen column density was allowed 
to vary freely (~1 x 10^^ cm~^), such a soft component still has ~96% of its X-ray fiux absorbed. 
If the hydrogen column density decreased dramatically — to Nh of a few times 10^^ cm~^ — the 
accretion shock emission could dominate the observed fiux. This seems a plausible scenario to 
explain the softer plasma detected in late 2004 March, especially if the large intervening hydrogen 
column density inferred at other observing epochs was due mostly to the accretion streams. 

An alternative explanation is that we observed V1647 Ori when it was in the midst of a large 
accretion episode that pushe d the star-disk boundary inward to the point where the accretion 



became non-magnetospheric (iHartmannl 1 1 9981 ) . effectively reducing the amount of hard X-ray fiux 



produced by magnetic reconnection events in the accretion stream. 

During the 2008-2009 epoch, when the X-ray luminosity of V1647 Ori increased, the spectrum 
hardened and the emitting plasma increased in temperature (Fig. [71 right panel); also, when the 
overall X-ray luminosity decreased, the X-ray spectrum softened and the X-ray generating plasma 
cooled. These correlations are also seen in the Chandra observations following the 2003 eruption 
(Fig. [Tj left panel). However, between 2008 September and 2008 November, the spectrum appeared 
to harden slightly as the X-ray luminosity decreased slightly. If we are observing X-rays generated 
predominantly by the ~1 keV plasma in the accretion footprint, then a decline in accretion should 
result in a decline in total X-ray luminosity, particularly in the lower energy fiux, resulting in 
a hardening of the spectrum. Therefore, the hardening of the spectrum and decrease in X-ray 
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luminosity observed between 2008 September and November might be a time when CXO was able 
to see plasma closer to the protostellar surface. In fact, this epoch showed a low temperature 
(~2 keV) plasma, consistent with what could be a mixture of hotter (4-6 keV) and cooler (1 
keV) plasmas. These data may offer evidence that the X-ray flux includes an emission component 
from the cooler plasma in the footprint of the accretion funnel located in the stellar photosphere. 
Because the cooler plasma suffers greater extinction and because CXO is less sensitive to softer 
X-ray emission than XMM-Newton, CXO would likely only detect the cooler plasma when the 
viewing geometry is favorable. 



4.1.1. Intervening Hydrogen Column Density 



Our modeling work for the X-ray observations obtained in 2004-2005 and 2008-2009, when 
the mean X-ray count rates were greater than 1 c ts ks~^ yield a best fit value for Nh of about 4.1 x 
10^^ c m~^, consistent wi th the results derived bv lCrosso et al.l ( 20051 ) and a visual extinction of Ay 
~ 20 (|Vuong et al.l 120031 ). We were unable to fit Nfj in our mod eling work for the low count-rate 
observations during the quiescent period in 2005-2006; however, lAspin. Beck. &: ReipurthI (120081 ) 
obtained a best fit value for Ay of 19 it 2, based on their optical, near-infrared, and mid-infrared 
observations obtained in February 2007, which w as also during the quiescent perio d (based on the 
optical and near- infrared photometry reported by lAspin. Beck. &: ReipurthI (120081 )). The observed 
value of Ay for February 2007 lends strong support to our use of A^^^- = 4.1 x 10^^ cm~^ for our 
modeling of the quiescent epoch observations. In addition, together these data suggest that Nh 
and, by implication. Ay, remained essentially unchanged as inferred from the X-ray observations, 
whether V1647 Ori was in the quiescent or elevated X-ray state. 

On t he other hand, a s see n at longer wavelengths, the extinction toward V1647 Ori has 
changed. Abraham et al. (2004) derive A y = 13 from near- i nfrared data obtained in 1998 by 
2MASS . During the outbu r st in 2004-2005, iBriceiio et all (I2OO4I ) found Ay = 8-10 on 2004 Febru- 
ary 18, iReipurth k Aspinl ^20oi ] found Aj = 1 26, Ah = 0.81, ari d Ak = 0.5 o n 2004 Feburary 
18, all of which are consistent with Ay of ~5 ([Becklin et al.lll978l ). IVacca et al.l (120041) obtaine d 
Ay = 11 from measurements of the 3.1 fim water band on 2004 March 9, and lOjha et al.l (j2006l ). 



who made optical and near-i nfrared observatio ns from 2004 into very late 2005, reported Ay of 
~5 during the 2004 outburst. lOjha et al.l ()2006l ) then re port that Ay increased to ~10 by the end 
of 2005, when V1647 Ori had dramatically faded, and lAspin. Beck, fc ReipurthI (120081 ) reported 
Ay was as high as 19 by early 2007. Clearly, the extinction, as measured via optical and near- 
infrared measurements, changed first from quiescence to outburst and then from outburst back to 
quiescence. 

These apparently discordant results have a straightforward and consistent explanation in the 
context of an accretion episode. In X-rays, we are essentially detecting V1647 Ori along a direct line 
of sight to the stellar photosphere. Our results therefore indicate that the absorption along this di- 
rect line of sight — which likely includes at least part of a thick circumstellar disk that is tilted about 
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30 degrees from edge-on (|Acosta-Pulido et alj 120071 ) — does not change significantly as a function 
of time, despite the evident changes in X-ray luminosity. The optical and near-infrared photons we 
observe, however, emerge from the near-photosphere environment of V1647 Ori along two paths. 
One path, along our direct line of sight to the photosphere, produces heavily reddened and extincted 
light. The second path takes photons nearly perpendicular to o ur line of sight, through a n evacu- 
ated polar cavity, where they then scatter into our line of sight (jAcosta-Pulido et al.l 120071 ). These 
photons are bluer and much less heavily extincted than the line-of-sight photons. When V1647 Ori 
is in the quiescent state (1998, late 2005-2007), we see a faint, reddened, heavily extincted source 
because the contribution from scattered light is minimal. During the outburst state (2004-2005; 
2008-2009), we see a brighter, bluer source because the contribution to the continuum of scattered 
photons is large. 



4.2. Possibility of a Second Plasma 



Though the plasma temperature of V1647 Ori strongly correlates with the X-ray luminosity 
and hardness ratio (Fig. [7j), it is unclear whether we are observing a single-component plasma 
that increases or decreases in temperature and thus causes the observed changes in X-ray lumi- 
nosity or if a second, lower-temperature plasma is also present and whose contribution to the 
total spectrum is o verwhelmed by the hotter tempera ture plasma during outbursts. Other sources, 
such as V710 Tau (jShukla et al.ll2008l ) and EX Lupi (jGrosso et al.ll20ld ) have spectra that can be 
model ed as two-teinperatu re plasmas, with one component fading as the star returns to quiescent 
levels. iGrosso et al.l (|2005l ) was able to model the XMM-Newton observation of V1647 Ori in 2004 
April with a single-temperature plasma but also found a better fit using a two-component model. 
However, this observation of V1647 Ori was at least twice as long as most of the Chandra ob- 
servations from 2004-2009. In addition, XMM-Newton is more sensitive to lower-energy photons 
than CXO/ACIS. Therefore, it may be the case that V1647 Ori does have two plasmas that con- 
tribute to the X-ray spectra but that ACIS was unable to consistently detect the lower-temperature 
component due to the inherit limitations of its design, the shorter observation times of the CXO 
observations, and the high degree of softer X-ray absorption by the intervening hydrogen column. 

In order to test whether CXO would be able to detect two distinct plasma temperatures in 5 ks 
and 20 ks observations, we simulated a two-component pl asma (using fakeit) with a normalization 
ratio and plasma temperatures identical to those found by lGrosso et al.l ()2005l ) and suffering extinc- 
tion by the same hydrogen column density. The simulated spectra were then convolved with the 
responses of the front-illuminated ACIS-I3 and back-illuminated ACIS-S3 CCDs. The simulations 
incorporated appropriately-scaled normalization parameters so that the resultant count rates were 
comparable to those of the 2005 and 2008-2009 CXO observations. Each of the simulated spectra 
was first fit with a single-component model and then with a two-component model. All of the 
20 ks observations were found to be fit better using two-component models with each fit yielding 
values for plasma temperatures and normalization that were close to the original simulation param- 
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eters. F-tests also suggested a slight improvement in t he model fits i f the 75%/25% normalization 
ratio between the low/high components, as found by iGrosso et alj (j2005l ). was forced. When fit 
with a two-component model, the simulated 5 ks observation of the more sensitive ACIS-S3 CCD 
yielded a best fit (using the 75%/25% normalization ratio) that converged to a model with two 
nearly-identical temperature plasmas, i.e., to a single-temperature plasma, and thus does not lead 
to a better fit. Thus, it appear s that CXO s h ould have the sensitivity to detect a two-component 
plasma such as that found by iGrosso et al.l (120051 ) if the observations are long enough and the 
lower-temperature plasma flux is strong enough compared to the higher-temperature plasma flux. 

The spectral models for each of the five Chandra observations obtained in 2008-2009 consis- 
tently converge to fits with a single-temperature plasma that, to within the errors, is ~4 keV. As 
seen with the correlation between X-ray luminosity and the mean hardness ratio, median photon 
energy, and mean count rate, we also find that when the X-ray luminosity increases or decreases, so 
does the plasma temperature. This single-compo nent model is consis tent with the single-component 
model of the 2005 April 11 XMM-Newton data (jCrosso et allbood ). 



4.3. Similarity of the 2003 and 2008 Eruptions 

Finally, the second (2008) eruption of V1647 Ori has very similar spectral characteristics to 
that of the first (2003-2005) eruption. From Tables [H El andU the five 2008-2009 observations of 
V1647 Ori, which span a period of approximately seven months, show 



• mean hardness ratios that are consistently as hard (to within the errors) as the hardest value 
of any observation following the previous eruption; 

• median photon energies at levels that are very similar to the greatest median photon energy 
of any observation from the previous eruption; 

• mean count rates that are 20-50 times higher than those observed during X-ray quiescence 
and that are typically 1-2 times the highest mean count rate of the 2003-2006 observations; 

• plasma emission measures that are usually 1-2 times greater than any modeled from the 
2003-2005 observations; 

• and X-ray luminosities that are consistently greater than 1.5 x 10^'^ ergs s~^, 1-2 times the 
X-ray luminosity of VI 647 Ori during any of the 2003-2005 CXO observations. 



Given that the spectral characteristics of VI 647 Ori are so similar in observations in which the 
X-ray fiux is elevated above the X-ray quiescent level, it appears that the same X-ray generation 
mechanism was at work during both eruptions and that the plasma characteristics were very similar 
during the two eruption epochs. The derived X-ray luminosities suggest that the second eruption 
was more energetic than the first, and given that the modeled emission measures are greater for 
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the observations following the second eruption, it is reasonable to conclude that we observed the 
same phenomenon in both eruptions but that a larger mass of X-ray emitting plasma was active 
during the outburst that began in 2008. 

5. Summary 

Our X-ray monitoring demonstrates that the two optical/ncar-IR outbursts undergone by the 
enigmatic V1647 Ori in less than a decade were accompanied by strong X-ray outbursts. During 
these two outbursts (2003-2005 and 2008-present), we see that the X-ray flux rose to peak lumi- 
nosity over a span of a few weeks and then remained elevated for approximately two years during 
the first eruption and for at least one year during the second eruption. Given that there is very 
strong evidence that the outbursts observed in the optical and near-infrared regimes are driven by 
accretion, we conclude that the correlated outbursts in X-rays are also driven by accretion. 

We find that all of the CXO spectra of V1647 Ori are best modeled with a single, moderate- 
temperature (2-6 keV) plasma. In almost all cases, the plasma temperature that emerges from 
models of the CXO spectra is too high to be generated via accretion hotspots on the stellar pho- 
tosphere but is reasonable for a plasma generated via magnetic reconnection events. However, 
the X-ray-emitting plasma could also be located in a strongly enhanced stellar corona, or at the 
inner edge of the circumstellar disk. Given that accretion is ongoing, lower-temperature plasma 
generated by shocks at the accretion footprint is very likely present; however, during these CXO 
observations, any such soft component contributed much less flux than the moderate-temperature 
plasma and so usually did not leave a distinct signature in the X-ray spectra. With the elevated 
hard X-ray flux lasting the duration of the 2008-2009 epochs, we conclude that the X-ray flux is 
not the result of typical coronal flares generated via reconnection events. We believe that since 
the optical/near-infrared flux remains elevated throughout this observation epoch, we arc instead 
observing the X-rays generated from reconnection events in the accretion stream, with the softer 
X-ray flux possibly being generated by accretion hotspots at the stellar photosphere. 

We find no significant change in X-ray absorbing column, indicating that varying optical/IR 
color mcasTircmcnts, which have previously been interpreted as evidence for variable reddening 
toward V1647 Ori, may instead be indicative of varying contributions from scattered vs. direct 
photospheric emission from the YSO. Two of the spectra obtained during the most recent (2008) 
outburst appear to show the 6.4 keV neutral iron feature indicating fluorescence from cold (pre- 
sumably circumstellar disk) gas surrounding V1647 Ori. 

With V1647 Ori being observed intensely during both outbursts at X-ray, optical, and infrared 
wavelengths, this objects stands as one of the best characterized systems that exhibits such a close 
correspondence between X-ray output and accretion rate. As a result of intense monitoring at 

X-ray, optical, and infrared wavelengths during two successive accretion-driven outbursts, V1647 
Ori stands as the best characterized YSO in terms of the correspondence between X-ray output 
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and accretion rate. We have shown, furthermore, that this YSO exhibited strikingly similar X-ray 

behavior and spectral properties during its recent accretion bursts. These results underscore the 
need for X-ray monitoring of additional eruptive YSOs, so as to evaluate whether the remarkable 
consistency of V1647 Ori is the exception or the norm. 

We thank Nuria Calvet for providing early access to the data in CXO ObsIDs 10763 and 
8585. This research was supported via awards numbers GO8-9016X and GO9-0006X to Vanderbilt 
University issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian 
Astrophysical Observatory for and on behalf of NASA under contract NAS8-03060. 
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Table 1. 2002-2009 ACIS Observations. 



ObsID 


Observation 


JD 


ACIS 


Exposure 


Net Counts 


Mean Count 


Median Photon 


Mean Hardness^ 


Hardness Ratio'' 




Date 




Chip 


(ks) 


(0.5-8.0 keV) 


Rate (ks-l) 


Energy (keV) 


Ratio 


of Total Counts 


2539'= 


2002 Nov 14 


2452593 


S2 


62.8 


17 


0.3 ± 0.1 


2.0 ± 0.3 


-0.34 ± 0.31 


-0.28 ± 0.28 


5307<^ 


2004 Mar 07 


2453072 


S3 


5.5 


64 


11.7 ± 1.5 


3.6 ± 0.3 


0.46 ± 0.15 


0.45 ± 0.14 


5308<^ 


2004 Mar 22 


2453087 


S3 


4.9 


12 


2.5 ± 0.8 


2.0 ± 0.5 


-0.45 ± 0.33 


-0.46 ± 0.33 


5382 


2005 Apr 11 


2453472 


13 


18.2 


86 


4.8 ± 0.5 


3.5 ± 0.1 


0.51 ± 0.19 


0.46 ± 0.12 


5383 


2005 Aug 27 


2453610 


13 


19.9 


18 


0.9 ± 0.2 


3.0 ± 0.8 


0.01 ± 0.29 


0.34 ± 0.25 


5384 


2005 Dec 09 


2453714 


13 


19.7 


2 


0.1 ± 0.1 


2.2 ± 1.1 


0.0 ± 0.20 


-0.07 ± 0.92 


6413 


2005 Dec 14 


2453719 


13 


18.1 


4 


0.2 ± 0.1 


1.7 ± 0.3 


-0.36 ± 0.26 


-1.0 ± 0.75 


6414 


2006 May 01 


2453857 


13 


21.6 


3 


0.1 ± 0.1 


1.3 ± 0.3 


0.09 ± 0.23 


0.07 ± 0.92 


6415 


2006 Aug 07 


2453955 


13 


20.5 


4 


0.2 ± 0.1 


2.4 ± 0.4 


-0.27 ± 0.21 


-1.0 ± 0.73 


9915 


2008 Sep 18 


2454728 


13 


19.9 


466 


23.6 ± 1.1 


3.5 ± 0.1 


0.36 ± 0.07 


0.37 ± 0.05 


10763'= 


2008 Nov 27 


2454798 


12 


19.7 


217 


11.0 ± 0.7 


3.9 ± 0.1 


0.65 ± 0.09 


0.64 ± 0.08 


8585= 


2008 Nov 28 


2454799 


12 


28.5 


160 


5.8 ± 0.4 


3.4 ± 0.1 


0.40 ± 0.10 


0.43 ± 0.09 


9916 


2009 Jan 23 


2454855 


13 


18.4 


245 


13.6 ± 0.9 


3.7 ± 0.2 


0.41 ± 0.08 


0.43 ± 0.07 


9917 


2009 Apr 21 


2454943 


13 


29.8 


260 


8.8 ± 0.5 


3.5 ± 0.2 


0.28 ± 0.08 


0.37 ± 0.07 



^Average of hardness ratios computed from 2 ks light curve data bins (ObsID 2539 uses 10 ks light curve data bins). 
''Hardness ratio computed using the total numbers of hard and soft X-ray photons from the entire observation. 

'=V1647 Ori was not the target of the observation, and given values have been adjusted to compensate for the 4' off-axis position of 
V1647 Ori. Exposure maps indicated net counts and mean count rates for ObsID 2539 should be increased by 8% while those for ObsIDs 
10763 and 8585 required a 3% increase. 

'^The back-illuminated ACIS-S3 CCD is more sensitive to X-rays from plasma in the temperature regime characteristic of V1647 Ori 
than are the front-illuminated ACIS CCDs. Values displayed for net counts and mean count rates, and their associated errors, for ObsIDs 
5307 and 5308 have been adjusted downward by 10%, based on spectral simulations in order to make net counts and mean count rates 
directly comparable between all chips. 

Note. — All errors are la. The net counts for each observation are the total number of counts w ithin the 0.5—8.0 keV range. Median 
photon energy uncertainties were calculated via the half-sample method used in lKastner et al" lHooi). Mean count rates were determined 
by dividing the net counts by exposure times. Uncertainties for mean count rates and hardness ratios of total counts follow Poisson 
statistics. The uncertainty for the hardness ratio of the total counts for ObsID 5384 could not be calculated due to the even distribution 
of the very low number of counts 
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Table 2. Comparison of Simulated Front- and Back-Illuminated Chip Exposures 



Plasma 


S3/I3 


S3 


13 


HR 


Temperature (kev) 


Count Ratio 


HR 


HR 


Difference 


2 


1.21 


0.02 


0.12 


0.10 


5 


1.14 


0.30 


0.38 


0.08 


7 


1.11 


0.40 


0.48 


0.08 



Note. — Net count and hardness ratio (of total counts) compar- 
isons for simulated spectra convolved ■with instrument responses from 
ObsIDs 5307 and 9915, which used the back-illuminated S3 CCD and 
front-illuminated 13 CCD, respectively. Each of the 1-ks spectra were 
simulated using the fakeit command in XSPEC with the intervening 
hydrogen column density and chemical abundance fixed at 4.0 x lO^^-^ 
cm~^ and 0.8 solar, respectively. The front-illuminated chip detects 
roughly 85% of the net counts of the back-illuminated chip, with a 
greater drop in sensitivity for soft versus hard X-rays. 



Table 3. Model Fits for 2008-2009 Chandra Observations. 



ObsID 


Observation 
Date 


Reduced 


Degrees of 
Freedom 


Nff 

(X 1022 cm-2) 


kTx 
(keV) 


EM 
(x 10^3 cm^ 


Observed Fx 
'-'') (x 10-13 ergs cm-2 ^ 


Observed Lx 
3-1) (x 1030 ergs s-i) 


9915 

10763 

8585 

9916 

9917 


2008 Sep 18 
2008 Nov 27 

2008 Nov 28 

2009 Jan 23 
2009 Apr 21 


0.85 
1.24 
0.71 
0.79 
0.88 


77 

39 
26 
42 
44 


4 1+0-9 
^•^ -0.7 

4.1 (Fixed) 

5 9+2-'' 
4 1+1-1 

-0.8 

4 1+1-3 

^•-^ -0.9 


-1.0 
4.0 (Fixed) 

3.41- 


15.0 1«-3 

6-1 ^si 

8.21^:^ 

6.8i?:« 

5 9+''-i 


^-^-1.1 

1 a+0.3 

10+0-1 

^-"-0.7 

2 6+°-2 

1 5+0-2 


8.11°:? 

1.8l- 

5.oi«:3 
2-9 1?:^ 



Note. — Uncertainties given for hydrogen column density, plasma temperature, and chemical abundance correspond to the 90% confidence 
intervals, whereas the observed X-ray luminosities are given with their corresponding 68% confidence intervals. Chemical abundance was 
fixed at 0.8 solar. For ObsID 10763, in order to constrain the observed X-ray flux and luminosity, the hydrogen column density and plasma 
temperature were fixed at the given values. Emission measures and luminosities assume a distance of 400 pc to V1647 Ori. 
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Table 4. Model Fits for 2002-2006 Chandra Observations. 



ObsID 


Observation 


Reduced 


Degrees of 




EM 


Observed Fx 


Observed hx 




Date 


Statistic 


Freedom 


(keV) 


(x 10^3 cm-3) 


(x 10~^^ ergs cm~^ s~ 


■1) (x 10^0 ergs s-l) 


2539 


2002 Nov 14 


1.37 


13 


09+^'^ 

"■■^ -0.5 


06+^'^ 


n. 03+0.01 


n5+° °l 


5307 


2004 Mar 7 


1.19=^ 


10 


4.0 (F) 


5.4 ± 1.3 


1.6 ± 0.5 


3.1 ± 0.9 


5308 


2004 Mar 22 


1.45 


10 


1 1+1-° 

^■^ -0.5 
3 4+3.0 


3 g+6.9 

3-1 ^1;^ 


0.2 ± O.l'' 


0.4 ± 0.1 


5382 


2005 Apr 11 


0.98=^ 


14 




1 5+°-2 


5383 


2005 Aug 27 


1.99 


16 


o 9+10.1 
■^■^ -1.7 


0.8+»;« 


0.2 


4+0.1 
"■^-0.3 


5384 


2005 Dec 9 


1.67 


1 


0.86 (F) 


Q 2+0-6 


0.01 ± O.Ol'' 


0.01 ± 0.01 


6413 


2005 Dec 14 


0.87 


3 


0.86 (F) 


9+°* 


0.02 ± O.Ol'' 


0.05 ± 0.03 


6414 


2006 May 1 


2.83 


2 


0.86 (F) 


3+0.6 


0.01 ± O.Ol'' 


0.02 ± 0.02 


6415 


2006 Aug 7 


1.45 


3 


0.86 (F) 


"•^ -0.3 


0.01 ± O.Ol'^ 


0.01 ± 0.01 



^Value is the reduced-x^ value. 

''Observed X-ray flux errors and corresponding X-ray luminosity errors were derived by multiplying the mean count rate 
errors from Table [T] by an energy conversion factor (EOF). 

Note. — Uncertainties given for hydrogen column density and plasma temperature correspond to the 90% confidence 
intervals, whereas the observed X-ray luminosities are given with their corresponding 68% confidence intervals. All models use 
a fixed hydrogen column density of N^f = 4.1 X 10^^ cm~^ and a chemical abundance fixed at 0.8 solar. For ObsID 5307, the 
plasma temperature was fixed ("F") at a typical post-outburst plasma temperature of kTx = 4.0 keV in order to constrain the 
remaining model parameters. For ObsIDs 5384, 6413, 6414, and 6415, the very low number of counts did not permit spectral 
fi tting with fr eely-varying plasma temperature, so the plasma temperature was fixed at a pre-outburst value as discussed in 
iKastner et al. j mo4 ). Observed fluxes and emission measures were derived from the resulting model fits, an energy conversion 
factor (EOF) was calculated for each observation, and the mean count rate uncertainties were multiplied with the EOF to 
better constrain the source fluxes displayed. All models, unless otherwise noted, use the Cash statistic to determine goodness of 
fit, and the reduced statistic is the Cash-statistic or value divided by the number of degrees of freedom. Emission measures 
and observed luminosities assume a distance of 400 pc to V1647 Ori. 
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Fig. 1. — : Background-subtracted X-ray light curves for 2002-2006 epoch observations. Time bins 
are 2 ks for each observation (except 2002 November 14, which uses 10 ks bins due to the very low 
count rate) and contain counts (0.5-8.0 keV) per total time associated with each bin, not the total 
time associated with the observation. Observations for 2005 December 09, 2006 May 01, and 2006 
August 07 do not detect V1647 Ori, so their light curves are not presented. Uncertainties in mean 
hardness ratios and count rates are la. 
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Fig. 2. — : X-ray light curves and hardness ratio of V1647 Ori during observation 9915. The top 
panel light curve covers the 0.5-8.0 keV energy range, the second panel down spans the lower 
energy range of 0.5-2.8 keV, and the third panel down covers the higher 2.8-8.0 keV energies. 
The bottom panel displays the hardness ratio for each of the bins. Counts were binned into 2 ks 
bins. Uncertainties in mean hardness ratios and count rates are Icr. Apparent "X-ray bright" time 
intervals are indicated by the shaded regions in the figure (4-8 ks and 16-20 ks) and in Figures 
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Fig. 3. — : X-ray light curves and hardness ratio of V1647 Ori during observation 10763. See 
Figure [2] for description of panels. Counts were binned into 2 ks bins. Uncertainties in mean 
hardness ratios and count rates are la. V1647 Ori had a slightly elevated X-ray flux, seen best in 
the hard X-ray band, from 8-10 ks. 
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Fig. 4. — : X-ray light curves and hardness ratio of V1647 Ori during observation 8585. Counts 
were binned into 2 ks bins. Uncertainties in mean hardness ratios and count rates are la. The 
hard and broad-band X-ray Ught curves show that the hard X-ray flux from V1647 Ori increased 
sharply at 20 ks into the observation and remained elevated for ~10 ks. During this period of 
elevated X-ray flux, the soft-band flux remained more or less constant. 
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Fig. 5. — : X-ray light curves and hardness ratio of V1647 Ori during observation 9916. Counts 
were binned into 2 ks bins. Uncertainties in mean hardness ratios and count rates are la. The 
X-ray hght curve shows that V1647 Ori began increasing in hard-band and, possibly, soft-band flux 
about 8 ks into the observation, and the flux remained elevated for ~12 ks. 
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Fig. 6. — : X-ray light curves and hardness ratio of V1647 Ori during observation 9917. Counts 
were binned into 2 ks bins. Uncertainties in mean hardness ratios and count rates are la. V1647 
Ori appeared to be in an elevated hard X-ray flux state at the onset of the observation, with the 
X-ray count rate remaining elevated for ~8 ks. 
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2002-2006 V1 647 Ori Chandra Observations 2008-2009 V1 647 Ori Chandra Observations 
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Fig. 7. — : Time series of observed X-ray luminosity (top panels), mean hardness ratio (middle 
panels), and plasma temperature (bottom panels) for V1647 Ori. Crosses represent data obtained 
with ACTS front-illuminated CCDs, and squares represent data obtained with the ACIS back- 
illuminated S3 CCD. Plotted uncertainties for the hardness ratios and luminosities represent the 
68% confidence interval (la) and plotted uncertainties for the plasma temperatures represent the 
90% confidence interval (1.6cr). The modeled plasma temperatures for the 2005 December through 
2006 May observations are not shown because this parameter was not well constrained by spectral 
fitting. 
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Fig. 8. — : Two-parameter plot for all CXO observations of V1647 Ori showing the correlation 
between the changes in mean hardness ratios of individual observations and the changes in observed 
X-ray luminosities. Crosses represent values corresponding to data obtained with ACIS front- 
illuminated CCDs, squares represent values involving data obtained with ACIS front-illuminated 
CCDs and the back-illuminated S3 CCD, and the triangle represents a value that used back- 
illuminated S3 CCD data only. The correlation coefficient for the changes in mean hardness ratio 
and changes in observed X-ray luminosity is 0.44. 
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Fig. 9. — : Best-fit XSPEC models of spectra obtained from 2008 September to 2009 April, with 
observations in chronological order from top to bottom. The top frame in each panel displays the 
data (binned to five-count-minimum bins) in black overlaid with the model in blue. 
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Fig. 10. — : Spectrum of the 2009 April 21 observation (black) overlaid with the single-component 
plasma model (blue) and the model with the additional Gaussian component (red) centered at 6.4 
keV (equivalent width of 510 eV) added to account for neutral iron emission. The entire modeled 
wavelength range (0.5-8.0 keV) is shown in the top panel, and the bottom panel displays the 5.0-8.0 
keV energy range to more easily show the iron emission and the Gaussian component. Spectral 
data points employ a three-count-minimum bin size. 
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V1647 Ori X-ray, Optical, and Infrared Variability 



2003 2004 2005 2006 2007 2008 2009 2010 2011 




2500 3000 3500 4000 4500 5000 5500 

Julian Date - 2450000 



Fig. 11. — : Near-infrared and X-ray light curves of V1647 Ori. Chandra X-ray data are 
shown as crosses; the Ic-band data are shown as s quares, and the H-ban d data are shown as 
circles. I-band Goddard Robotic Telescope (GRT) ( Sakamoto et al.ll201ll ) data (Michael Rich- 
mond, personal communication) were obtained in 2009. For the published I- and H-band 
data, e rror bar size is on th e order of the dat a point size (da t a are fr om: 



McGehee et al 



Oiha et al.1 (12009) 



Fedele et al 



(2007), 



Briceho et al 



(12004). iReipurth fc Aspinl (120041) ISemkovl (120041). lOiha et al 



Semkov (2006), Venkat fc Anandarao (2006), 



Acosta-Pulido et al 



Aspin Beck, k Reipurthl (|2008l ). lOiha et al.l (12003 ). lAspin et al.l (|2009l ). and 



2005 ). 



( 20071 ). 



Venkat &: Anandaraol (|201ll )). The uncertainties for the GRT data and CXO data are la. Calendar 



year is indicated along the top horizontal axis. 
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Fig. 12. — : Two-parameter plot for all CXO observations of V1647 Ori showing the correlation 
between the observed Ic-band near-infrared and X-ray luminosities. Crosses represent data obtained 
with ACIS front-illuminated CCDs, and squares represent data obtained with the ACIS back- 
illuminated S3 CCD. Errors for Ic-band luminosities that were interpolated at X-ray observation 
dates were calculated by averaging the errors of the Ic-band luminosities of the five observations 
made nearest to the interpolation date. Circles represent four observations (ObsIDs 9915, 10763, 
8585, and 9917) in which the interpolated Ic-band luminosity of ObsID 9916 was assumed to be the 
correlated Ic-band luminosity. The correlation coefficient for the Ic-band and X-ray luminosities 
without the circle values is 0.65, while the correlation coefficient using all values is 0.67. 



